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Abstract 

Background: Lupulin glands of hop produce a specific metabolome including hop bitter acids valuable for the brewing 
process and prenylflavonoids with promising health-beneficial activities. The detailed analysis of the transcription factor 
(JF)-nnediated regulation of the oligofamily of one of the key enzymes, i.e., chalcone synthase CHS_H1 that efficiently 
catalyzes the production of naringenin chalcone, a direct precursor of prenylflavonoids in hop, constitutes an important 
part of the dissection of the biosynthetic pathways leading to the accumulation of these compounds. 

Results: Homologues of flavonoid-regulating TFs HIMyb2 (M2), HlbHLH2 (B2) and HIWDR1 (W1) from hop were 
cloned using a lupulin gland-specific cDNA library from the hop variety Osvald's 72. Using a "combinatorial" 
transient GUS expression system it was shown that these unique lupulin-gland-associated TFs significantly activated 
the promoter (P) of chs_H] in ternary combinations of B2, W1 and either M2 or the previously characterized 
/-//Myb3 (M3). The promoter activation was strongly dependent on the Myb-P binding box TCCTACC having a core 
sequence CCWACC positioned on its 5' end region and it seems that the complexity of the promoter plays an 
important role. M2B2W1 -mediated activation significantly exceeded the strength of expression of native chs_V\] 
gene driven by the 35S promoter of CaMV, while M3B2W1 resulted in 30% of the 35S:c/7S_H1 expression level, as 
quantified by real-time PCR. Another newly cloned hop TF, /-//Myb7, containing a transcriptional repressor-like motif 
pdLNLD/ELxiG/S (PDLNLELRIS), was identified as an efficient inhibitor of chs_H] -activating TFs. Comparative 
analyses of hop and A tholiono TFs revealed a complex activation of Pchs_H] and Pchs4 in combinatorial or 
independent manners. 

Conclusions: This study on the sequences and functions of various lupulin gland-specific transcription factors 
provides insight into the complex character of the regulation of the c/75_H1 gene that depends on variable 
activation by combinations of R2R3Myb, bHLH and WDR TF homologues and inhibition by a Myb repressor. 



Background 

Hop {Humulus lupulus L.) plants are mainly cultivated for 
the brewing industry, as a source of flavor-active secondary 
metabolites contained in the lupulin glands, i.e. glandular 
trichomes that develop in the hop female inflorescences 
(cones). In addition, hop has been known for a long time 
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in traditional medicine and recently several compounds in 
the lupulin metabolome, including hop bitter acids and 
prenylated flavonoids, have received particular attention in 
view of their highly interesting medicinal properties e.g., 
[1-3]. Xanthohumol (X), the principal prenylated chalcone 
in the lupulin glands is a fascinating cancer- chemopreven- 
tive compound exhibiting a broad spectrum of inhibition 
mechanisms at all stages of carcinogenesis [4]. Although X 
was shown to be poorly bioavailable [2,5-7], interesting 
anti-inflammatory in vivo results have been obtained in 
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specific target tissues such as the liver [8]. Another lupulin- 
derived prenylflavonoid, 8-prenylnaringenin (8-PN), is one 
of the most potent phytoestrogens known to date [9,10]. 
Together with its precursor isoxanthohumol which can be 
metabolized in the body by the gut microbiota to 8-PN 
[11], these compounds are indicated as the active ingredi- 
ents of hop extracts targeting relief of menopausal symp- 
toms [12,13]. 

The biosynthesis of naringenin chalcone in hop cones as 
a prenylflavonoid precursor is attributed to the gene encod- 
ing the CHS_H1 protein having a so-called "true" chalcone 
synthase (EC 2.3.1.74) activity, by which it efficiently cata- 
lyzes the production of naringenin chalcone by condensa- 
tion of three malonyl-CoA units and /?-coumaroyl-CoA 
[14,15]. Moreover, recombinant CHSJHl can utilize isova- 
leryl and isobutyryl CoA substrates, albeit at a low rate, 
and, therefore, could also be involved in the biosynthesis of 
hop bitter acids [15]. CHS_H1 is encoded by an oligofamily 
of genes having very specific expression in hop cones [16]. 
Other chalcone synthase-like enzymes associated with 
lupulin glands have been described. For instance, valero- 
phenone synthase (VPS) [17] is an enzyme which may be 
considered as a major component in the biosynthesis of 
hop bitter acids. The VPS-like homologue CHS4 also 
shows high expression in lupulin glands, however the 
enzyme does not catalyze the formation of naringenin chal- 
cone and its function is currently still unknown [15,18]. 

The complexity of the promoter elements of the chs_Hl 
genes suggests the involvement of several types of tran- 
scription factors (TFs), mainly from Myb, bHLH, and bZip 
families [16,19-21], in either independent or combinatorial 
pathways [22]. Independent or combinatorial activity of 
TFs in the regulation of flavonoid biosynthetic pathways 
has been described in several recent reviews [23-26]. 
Three TFs, i.e., R2R3Myb (M), bHLH (B), and WDR (W), 
exert combinatorial activation by formation of ternary 
complexes (MBW complexes) through proteimprotein 
interactions [26]. Such MBW complexes are highly orga- 
nized and each subunit fulfills a specific function such as 
binding to DNA, activation of expression of a target gene 
or stabilization of the transcription factor complex [23]. 
MBW complexes have been clearly identified in the flavo- 
noid biosynthetic pathway of Arabidopsis thaliana and 
Petunia hybrida, respectively, as a TT2/TT8/TTG1 com- 
plex (Transparent Testa 2/Transparent Testa 8/Transpar- 
ent Testa Glabra 1) driving coloration of the seed coat 
[27,28] and a triple combination complex AN2/AN1/ 
AN11 (Anthocyanin 2/1/11) regulating anthocyanin accu- 
mulation in the corolla; for reviews, see also [23,26,29]. 
The regulation of the flavonoid biosynthesis pathway by 
ternary complexes has also been shown in peas [30] and 
Lotus japonicus [31] and has been proposed for Perilla fru- 
tescens [32], Japanese morning glory [33], grapevine [34] 
and Pyrus [35]. 



Besides triple combinations, numerous examples from 
classic genetic or molecular genetic studies using ectopic 
TF expression analysis have documented or predicted 
enhanced flavonoid biosynthesis by binary complexes 
composed of two TFs from the R2R3Myb, bHLH, and 
WDR classes [34,36-40]. 

For the hop chs_Hl genes, direct promoter activation 
by the heterologous PAP1 TF from Arabidopsis thaliana 
and H/bZIPl and 7//bZIP2 TFs from hop has been 
demonstrated previously using a transient expression sys- 
tem [16,21]. Two hop R2R3Myb TFs, i.e. H/Mybl and 
///Myb3, have been characterized in our previous studies 
[19,20] and diverse biological effects caused by ///Myb3 
subvariants have been demonstrated in heterologous 
transgenotes [20]. Although these Myb TFs were sug- 
gested to be involved in the lupulin metabolome produc- 
tion based on sequence similarity and specific expression 
in hop cones [19,20], their influence on the activation of 
chs genes has not been investigated in detail. 

In the present work, we cloned novel lupulin-specific 
TFs and showed by a combinatorial transient expression 
assay that the 7//bHLH2 and MWDR1 TFs strongly acti- 
vate the chs _H1 genes in combinations with HlMyb2 and 
HlMyb3 TFs from hop, suggesting the formation of tern- 
ary complexes. This study confirms that the mode of 
action of the TFs strongly depends on the composition of 
the chs promoter and that HlMyb7 acts as a repressor of 
activating complexes. 

Results 

Cloning, comparative sequence and genomic analyses of 
R2R3Myb, bHLH and WDR transcription factors from hop 

During our previous work, we cloned two authentic 
R2R3Myb TFs having specific expression in cones of the 
Czech hop variety Osvald's 72, i.e., HIMybl [GenBank: 
AJ876882] [19] and MMyb3 [GenBank:AM501509] [20]. 
In the present work, we isolated and purified lupulin glands 
from Osvald's 72 according to Nagel et al. [41] and con- 
structed a lupulin gland- specific cDNA library. This cDNA 
library is more specific and more complex (3.5 x 10 6 pfu) 
than the previously constructed cDNA library from hop 
inflorescences and cones [19]. We used sequence motifs 
for two R2R3Myb TFs available in the GenBank database 
[no. AB292245 and AB292244], to amplify HlMyb2 [Gen- 
Bank: FN646081] and HlMyb7 [GenBank:FR873650], 
respectively, from the cDNA library. Both clones differed 
in several point mutations and, in addition, an 18 bp dele- 
tion was found for the clone HlMyb2 in comparison with 
the original database sequence [GenBank:AB292245]. 
The full-length clones of HlbHLH2 [GenBank:FR751553] 
were obtained by screening of hop EST sequences after 
pre-selections of 5' and 3' cDNA motifs in the TrichOME 
database (http://trichome.noble.org/trichomedb/). To 
get full length cDNA clones of H1WDR1 [GenBank: 
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FN689721] we firstly identified the 5' region from EST of 
the H. lupulus var. Phoenix. Full-length cDNA was 
obtained from a lupulin gland-specific library of the 
Osvald's 72 cv. using a combination of a 5primer and an 
oligo-dT anchor primer followed by nested PCR as 
described in Methods. The authenticity of all newly 
obtained cds sequences was verified by sequencing of the 
high-fidelity RT-PCR clones. 

At present, 7 full-length R2R3Myb factors from the 
Osvald's 72 cultivar are available in the sequence database. 
From these sequences, the two newly cloned ///Myb2 and 
///Myb7 showed high amino acid identity from 67 to 82% 
within the R2R3 domain to A. thaliana AtMybl2 [Gen- 
Bank: NM_130314] and PI TFs from Z. mays [GenBank: 
NM_001111873] involved in the flavonoid biosynthesis. A 
more wide comparison with various related R2R3 Myb TFs 
selected from BLAST analysis shows the positioning of 
HlMyb7 and ///Myb2 in different closely related clusters 
(Figure 1). The previously cloned ///Mybl [19] is clustering 
together with ///Myb7, however, this sequence is less 
related to ZmPl. HlMyb?> that was found earlier as metabo- 
lome regulator [20] is clustering together with C. sativa, V. 
vinifera, and G. hirsutum Mybs as described previously 
[20]. This TF is clustering together with ///Myb6, while two 
other hop Mybs 4 and 5 are quite unrelated and are not 
included in the analyzed dataset of Myb sequences (see 
Additional file 1). Within the presented R2R3Mybs, some 
of them were proven to be involved in anthocyanin meta- 
bolism or/and trichome development forming ternary 
BMW complexes such as AtTT2 [GenBank: NM_122946], 
AtMyb75 [GenBank: NM_104541], AtMyb23 [GenBank: 
NM_123397] and PiAN2 from P. integrifolia [GenBank: 
AF146703]. 

///Myb2 represents a typical R2R3Myb, with a pre- 
dicted mass of 29.7 kDa and a pi of 8.16. ///Myb7 is an 
R2R3Myb with a predicted protein mass of 29.6 kDa and 
pi of 8.19. These characteristics are close to other Myb 
TFs forming ternary complexes like TT2 from A. thali- 
ana (29.6 kDa, pi 8.8) or AN2 from P. hybrida (29.0 kDa, 
pi 7.23). In the same range of the molecular mass, but 
more acidic, are the subvariants of hop Myb3, i.e., 
s-///Myb3 (29.9 kDa, pi 6.4) and 1-///Myb3 (30.3 kDa, pi 
6.19). In the C-terminal region of ///Myb7, we identified 
a PDLNLELRIS motif which conforms to the consensus 
pdLNLD/ELxiG/S characteristic for plant Myb repressors 
[42] described in A. thaliana and Z. mays [43,44] and 
showing similarities to II AP2/ERF transcriptional repres- 
sors [45]. While the R3 repeat in 7//Myb2 perfectly 
matches the sub-group Illf bHLH interaction motif [DE] 
Lx2[RK]x3Lx6Lx3R [24,38,46], the HIMybl sequence 
matches this motif only partly (66% identity) and in 
///Myb3 this motif is absent (see Additional file 2). How- 
ever, the following consensus sequence. LllcLHphhGN- 
RWShIAth_(according to Bork consensus, http://coot. 



embl.de/Alignment//consensus.html) is characteristic for 
this interacting region, where most positions of putative 
interacting residues (underlined) are hydrophobic amino 
acids exposed on the surface. 

The cloned HWHLH2 cDNA sequence encoded for an 
acidic (pi 5.1) protein having a molecular weight of 77.1 
kDa. According to analysis of bHLH domains of selected 
plant bHLH TFs sequences from BLAST and from the 
dataset of Carretero-Paulet et al.[47], HlbWUM shows the 
highest similarity to bHLH AN1 from P. hybrida (90% 
amino acid identity within the bHLH domain). It is cluster- 
ing with several bHLH TFs that are involved in the regula- 
tion of the anthocyanin pathway, and belongs to wide 
cluster covering TFs with biological functions connected to 
flavonoid/anthocyanin metabolism, trichome, stomatal 
complex and flower development within functional classes 
1, 2, 5, and 10 as arranged by Carretero-Paulet et al. [47] 
(Figure 2A). Some of the related sequences were described 
to form MBW complexes (Figure 2A, more in detail 
described in Additional file 3). Comparisons of the bHLH 
basic regions revealed the presence of amino acid motif 
(H9, E13, R17) in the 7//bHLH2 domain characteristic for 
G-box binders [47], similarly to well characterized bHLH 
AN1 from P. hybrida and bHLH TT8 from A. thaliana 
[48] (Figure 2B). 

The cloned HIWDR1 cDNA encodes for a 38 kDa protein 
having a pi of 4.9. Phylogenetic analysis of this TF using 
sequences selected from BLAST revealed closely related 
WDR clusters within the group of dicotyledonous plants 
(Figure 3 A, for more details see Additional file 4). Most of 
these TFs were proven to regulate anthocyanin metabolism 
and some of them to form MBW complexes. 7//WDR1 
showed the highest similarity to recently reported 
sequences isolated from pear, apple, plum and raspberry 
(Figure 3A). Hop WDR1 shares 78.6% of identical amino 
acid with the AN11 protein from P. hybrida originally char- 
acterized as a plant WDR by de Vetten et al. [49] (Figure 
3B). More detailed comparisons of 7//WDR1 with P/zANll 
identified five polypeptide sequences that constitute the 
WD-repeats according to [49]. Only nine non-equivalent 
positions were identified in the 7//WDR1 repeats in com- 
parison to P/zANll (Figure 3B, sequences underlined). 

In order to assess the complexity of the genes encoding 
the newly cloned hop transcription factors ///Myb2, 
HlbHUYl and 7//WDR1, Southern blot analysis was car- 
ried out on genomic DNA isolated from the hop cultivar 
Osvald's 72 and digested with several restriction enzymes 
selected on the basis of the known cDNA sequences. 
Hybridization with probes created from cDNA fragments 
of ///Myb2, HlbttUYl and ///WDR1 genes showed single 
bands for all tested genes (Figure 4) when digested with 
EcoRl (HlbHLH2 and HIWDR1), Xhol (HlbHLH2 and 
HlMyb2), Xbal (HlMyb2 and HIWDR1), Pstl (HlbHLH2), 
Bglll (HlMyb2) and BamWl (HIWDR1) enzymes. This 
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Figure 1 Phylogenetic analysis for amino acid sequences of selected plant MYB TFs. The sequences of R2 and R3 domains of 40 selected 
plant MYB TFs were compared with PhyML v.3.0 and the tree was rooted along /\fMyb>075 using FigTree v.1.3.1. (see Methods). The lengths of 
interrupted branches are expressed in aLTR values corresponding to the scale. See additional file 1 for amino acid alignment, ID numbers and 
annotations of individual TFs. 



suggests that the cloned hop transcription factors 
HlMyb2, HlbHim and HIWDR1 are rather unique in 
the hop genome of Osvald's 72. The occurrence of two 
bands in the case of HWHLH2 cleaved with EcoRl (5.7 
and 0.6 kb) and HlMyb2 cleaved with Xbal (6.0, 4.0 and 
3.2 kb) is likely due to cleavage within some intron(s). 

Comparative expression level analyses of the 
transcription factors HIMyb 2/7, HlbHLH2 and HIWDR1 in 
hop tissues show specificity for lupulin glands 

Myb, bHLH and WDR genes from hop were amplified 
by the use of sequence motifs preselected from tri- 
chome databases and cloned from a cDNA library pre- 
pared from purified glandular trichomes. This 



corroborated to an initial specificity for lupulin glands. 
In order to verify the specificity of the expression of 
the cloned hop TFs HlMyb2, HlMyb7, HWHLH2, 
HIWDR1 we analyzed mRNA levels in various tissues 
using Real-Time quantitative PCR (RT qPCR). Glycer- 
aldehyde-3-phosphate dehydrogenase (GAPDH) mRNA 
was selected as a housekeeping gene for normalization 
based on previous experience and reports from others 
[41,50,51]. In order to avoid RNA contaminations from 
hop cone bracts, the lupulin glands of clone Osvald's 
72 were carefully separated and purified as recom- 
mended by Nagel et al. [41]. The GAPDH-normalized 
results show that HlMyb2 expression is highly specific 
for lupulin glands (Figure 5A) as the relative expression 
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Figure 2 Phylogenetic analysis for amino acid sequences of selected plant bHLH TFs (A) and Alignment analysis of the bHLH2 domain 

(B). The sequences of bHLH domain of 58 selected plant bHLH TFs were analyzed using PhyML v. 3.0. The polar tree cladogram presented in 
panel A was visualized using FigTree v.1.3.1. (see Methods). The lengths of branches are expressed in aLTR values corresponding to the scale. 
Alignment of the bHLH domains of hop H/bHLH2, petunia PhAN1 and A thaliana AtbHLH042 (TT8) is shown in panel B. The bHLH domain of 
H/bHLH2 was identified by similarity to the TTS domain as described by [48]. Amino acids within the basic region crucial for G-box DNA binders 
[47] are indicated by squares. See additional file 3 for amino acid alignment, ID numbers and annotations of individual TFs. 
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Figure 3 Phylogenetic analysis for amino acid sequences of selected plant WDR TFs (A) and comparison with the predicted H/WDR1 
sequence (B). Twenty five selected entries of plant WDR TFs (for the data set see additional file 4) were analyzed using PhyML v. 3.0 and the 
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Figure 4 Southern blot analysis of 5 ug genomic DNA isolated 
from the hop cv. Osvald's clone 72. The blot was hybridized 
with a 32 P-labeled cDNA fragment. (A) HA/VDR1 gene: BomH\ (Ba) 
and Xba\ (Xb) do not cut within the H/WD40 gene, EcoR\ (E) has one 
restriction site 146 bp from the 5'end of the cDNA. (B) H/bHLH2 
gene: EcoR\ (E), Xho\ (Xh) and Pst\ (?) do not cut within the H/bHLH2 
gene. (C) H/Myb2 gene: Xba\ (Xb) does not cut in the H/Myb2 gene, 
BglW (Bg) and Xho\ (Xh) have one restriction size, but only 69 resp. 
43 bp from the 3'end of the cDNA. The arrows show the position of 
the restriction fragment hybridizing with the cDNA probe fragment. 
The estimated size is indicated in kilobases (kb). 



levels in flowers and young cones only amounted to 
approximately 5% while in other tissues only traces of 
the PCR product were found. Almost identical results 
were obtained for HlbHLH2 (Figure 5C), where expression 
of only about 1% were detected in tissues other than puri- 
fied lupulin glands. A 4% relative expression level was 
detected in whole young cones. The expression of 
HIWDR1 was less specific for lupulin glands (Figure 5D). 
A highly significant relative expression level of this TF, 
reaching more than 20%, was observed in petioles, flowers 
and in whole young cones, where maturating lupulin 
glands comprise only a negligible fraction of total weight 
of the cone tissues. In addition, immature pollen was 
found to have HIWDR1 mRNA relative expression levels 
of about 10%. Furthermore, lower expression of HIWDR1 
was detected in somatic tissues such as young leaves and 
roots (Figure 5D). The lowest specificity for lupulin glands 
was observed for HlMyb7, where roots and whole young 
cones were found to contain a relative expression level of 
57% and 35%, respectively, whereas other tissues including 
leaves, flowers and immature pollen also showed high 
expression levels of over 15% (Figure 5B). 

Combinatorial transient expression analysis of selected 
lupulin gland TFs revealed a complex activation of 
pchs_H1 and cfo_H1 

According to our results, HlMyb2 and HWHLH2 exhibited 
high specific expression in lupulin glands similar to the 
previously investigated HlMyb3 [21]. Although the expres- 
sion of HIWDR1 and HlMyb7 showed a lower specificity, 
their relative expression levels were still most pronounced 
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Figure 5 Real Time RT PCR analysis of hop TFs. HIMyb2 (A) 7 
HIMyb7 (B), HlbHLH2 (C) and HIWDR1 (D) expression in lupulin 
glands and other tissues of H. lupulus Osvald's clone 72 as 
designated in the figure. GAPDH was used as house-keeping 
gene. 



in the glandular trichomes. Given the high sequence simi- 
larity of all these TFs to known regulators of the flavonoid 
pathway and their specific expression pattern, it is likely 
that these TFs are involved in the regulation of genes like 
c/z5_Hl, that are anticipated to be involved in the bio- 
synthesis of secondary metabolites in the lupulin glands 
[52]. 
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To assess the function of the cloned TFs we have used 
a'combinatorial transient expression'system that was 
developed previously [16,21]. In this assay the activation 
of chs_Hl promoter (Pc/zs_Hl) fused to the GUS refer- 
ence gene (Figure 6) is measured after co-infiltration of 
A. tumefaciens bearing Pc/zs_Hl and hop TFs vectors in 
various combinations into N. benihamiana leaves. Only 
negligible background signal resulting from activation of 
VchsJHl by internal leaf TFs was detected [16,21]. 

Pc/zs_Hl contains several Myb and bHLH-binding 
boxes as described previously [16] and analyzed further 
in this work (Figure 6A). According to the PLACE data- 
base the sequence TCCTACC (MybPZm box) containing 
the consensus of the Z. mays P factor binding box 
CCWACC [53] was identified at position -481 to -476. 
This potential TF binding box was of high interest in 
relation to the significant similarity of the newly cloned 
HlMyb2 and HlMyb7 to this ZmPl factor (Figure 1). At 
position -428 to -412 a Myb-like box was identified using 
the Genomatix database (http://www.genomatix.de) with 
a consensus for GA-Myb proteins [54]. For position -336 
to -328 the sequence CACCAAAC conforming to the 
Myb-P-related consensus sequence MACCWAMC was 
found using the PLACE database [55] and at position 
-287 to -271 a plant MYB IIG-type-like binding site was 
identified, i.e., TAAGGTAGTTGA (mismatch under- 
lined), which is similar to the Z. mays CI Myb-domain 
protein (Myb-C box) [56]. In addition, two typical bHLH 
G-binding boxes with the sequence CACGTG, that could 
be important for 7//bHLH2-mediated activation of 
?chs_Hl [57], were identified at positions -224 to -229 
and -313 to -318 (Figure 6A). 

The maximum values of transient GUS activity was 
reached at 66 to 92 h post-infiltration and declined there- 
after dependent on the TF or TFs combination co-infil- 
trated. For example, HlMyb3 variants and all complexes 
containing the HlMyb3 gene peaked at 4 days post infiltra- 
tion (dpi), whereas the HlMyb2 gene reached maximum 
activity 66-68 h dpi, except for the M2W1 binary complex. 

According to our results, a weak activation of Vchs _H1 
was observed after co-infiltration with ///Myb 2 (Figure 3). 
A very weak activation signal on the level of 40 pmol MU/ 
mg/min was found for ///Myb3 in accordance to our pre- 
vious study [21]. A similar very low activation of VchsJHl 
was found for 7//bHLH2, despite the presence of two 
bHLH G-boxes (Figure 6A), and no activation exceeding 
the background signal was observed when 7//WDR1 TF 
was applied by itself (not shown). However, dramatic 
increases of the activation of Pc/zs_Hl were observed both 
in Myb/bHLH binary and ternary combinations of the TFs 
under investigation (Figure 7A). The ternary combination 
led to 1721 ± 315 pmol MU/mg/min for M3B2W1 and 
5089 ± 107 pmol MU/mg/min for M2B2W1 (Figure 7A). 
The M2B2W1 combination significantly exceeded the 



activation caused by the strong constitutive 35S promoter 
of CaMV (1658 ± 102 pmol MU/mg/min). 

In order to assay the real activation of transcription of 
chs_Hl genes having very conserved promoter regions 
[16] we performed the co-infiltration of TFs with the 
"native" chs_Hl gene containing the natural intron 
sequence (Figure 6B) as described earlier [16]. The relative 
expression levels of chs_Hl mRNA were then measured 
by RT qPCR and 35S:c/zs_Hl was used as a reference (Fig- 
ure 7B). In correspondence to the previous results, the 
ternary combinations M3B2W1 and M2B2W1 were found 
to strongly activate the natural chs_Hl gene construct 
thereby showing good accordance between both assays. 

To enable the characterization of the binding sites 
responsible for the strong activation observed with both 
M3B2W1 and M2B2W1, a series of truncated variants of 
Pc/z5_Hl lacking the predicted TF-binding boxes were 
inserted (Figure 6A) in the GUS reference vector (Figure 
6B) and assayed in a transient expression system (Table 1). 
The majority of promoter activity (about 50%) was found 
to depend on the presence of the MybP box, although a 
significant activity (about 50%) was still observed in variant 
2 containing the Myb-like box, Myb-P like box, Myb-C 
box and two G-boxes (Figure 6 A, Table 1). By removal of 
the boxes, the promoter activity is gradually decreasing in 
a similar fashion for both ternary combinations to the low 
activity of construct 5 containing a single G-box. In order 
to confirm the high dependency of Pc/zs_Hl on the Myb-P 
box, we changed the sequence TCCTACC by a triple 
mutation to T AAC ACC (mutated nucleotides underlined) 
(Table 1). A decrease of about 50% was observed for both 
ternary combinations. Despite the critical role of the Myb- 
P box, a APc/z5_Hl variant containing the Myb-P box 
attached upstream of the TATA box (Figure 6A) was inac- 
tive, suggesting the dependence of VchsJM activation on 
the length and/or specific order in the promoter sequence. 

Substitution and complementation of activating 
complexes with selected H. Lupulus- and A Thaliana- 
derived myb TFs 

To investigate a possible functional complementation, 
related Myb TFs from H. lupulus and A. thaliana were 
compared using the transient expression system. In addi- 
tion to Pc/z5_Hl, also the promoter Pchs4 of the chalcone 
synthase-like gene from hop was analyzed. The gene for 
chs4 was cloned previously [GenBank: AJ430353] by Novak 
et al. [18]. The Pc/zs4 sequence is 629 bp in length and has 
a different architecture in comparison to Pc/zs_Hl, while 
still containing similar Myb and bHLH boxes. Using the 
PLACE database it was possible to identify a MybPZm box 
core CCWACC (sequence CCAACC) [53] at position -101 
to -95, a Myb P-like box corresponding to the consensus 
sequence MACCWAMC (sequence AACCTAAC) at posi- 
tion -92 to -85 [55] and another potential MyblAT box 
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Figure 6 Schematic drawing of the analyzed variants of the d7s_H1 promoter (A) and expression cassettes within the T-DNA parts of 
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Figure 7 Analysis of activation of the d7s_H1 promoter (A) and the whole d7s_H1 gene (B). The combinatorial transient expression assay 
by analysis of GUS activity was performed to analyze promoter activation (A); the relative levels of c/is_H1 mRNA were determined by RT qPCR 
and normalized against the elongation factor 1 housekeeping gene to analyze the activation of c/is_H1 gene by the TFs (B). Infiltrated A 
tumefaciens strains in panel A: C-control, infiltration of A tumefaciens LBA 4404 without plant vector; M2, H/Myb2; s-M3; s-H/Myb3; +W1, 
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with the consensus sequence WAACCA [59] at position 
-212 to -207 (sequence TAACCA). One bHLH G-binding 
box (sequence CACGTG) and one E-binding box 
(CACATG) [60] were found at positions -132 to -127 and 
-120 to -115, respectively. 

The relative activation levels for the individual TF 
combinations are shown in Table 2. It was observed 
that the Pc/zs_Hl promoter is most strongly responding 



to the ternary combinations, whereas the strongest 
response of Pc/zs4 was found for independent Myb action 
or binary combinations. Furthermore, A. thaliana-derived 
A£Myb23 {At5g40330) andA£Myb75 (AflPAPl, A tlg56650) 
clearly supplemented the M component in the ternary 
MBW combination, while A£Mybl2 did not cause any 
enhancement of activation in the MBW combination over 
the MB or MW combinations. A£Mybl2 also strongly 
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Table 1 Activation of modified variants of Pcfrs_H1 by 
hop TFs § 



Promoter variant* 


length (bp) 


Combinations 
M2B2W1 


M3B2W1 


P]chs H1 


500 


1 00.0*** 


1 00.0*** 


P Mut chs_H1** 


500 


56.3 ± 26.0 


48.0 ± 8.8 


P2c/7S_H1 


436 


52.2 ± 16.7 


40.7 ± 22.1 


P3c/is_H1 


346 


41.7 ± 18.8 


21.1 ± 2.0 


P4c/is_H1 


297 


12.2 ± 2.3 


12.2 ± 7.0 


P5c/7S_H1 


240 


3.5 ± 0.6 


2.1 ± 0.1 


P6c/is_H1 


217 


0.4 ± 0.1 


0.3 ± 0.1 


PAc/7S_H1 


271 


0.7 ± 0.2 


0.6 ± 0.1 


BGF-0* 


0 


« 0.1 


« 0.1 


P35S:BGF ## 


681 


29.2 ± 15.6 


87.5 ± 5.3 



§ Supplemental GUS histochemical analysis together with graph containing 
GUS activity values is presented in additional file 8. 

* for details see Figure 6A; ** Myb-P box was mutated; *** 100% was 5089 ± 
84 and 1721 ±315 pmol MU/mg/min for M2B2W1 and M3B2W1, respectively 

# infiltration with promoterless BGF-0 vector [58] (negative control, 
background); ## infiltration with p35S:BGF vector (35S control - see Methods); 
the average activity measured for p35S:GUS without co-infiltration with MBW 
TFs was 1658 ± 102 pmol MU/mg/min 

activated Pchs4 when applied alone or in the MW combi- 
nation. A£Myb75 and ///Myb2, the strongest ternary acti- 
vators of Pc/z5_Hl, showed low capability to activate 
Pchs4. Interestingly, a clear-cut difference was observed 
for the subvariants of HlMyb3 [20]. While /-///Myb 3 
shows a maximum of activation for Pc/zs4 in the MB com- 
bination, the maximum of activation of Pc/zs_Hl was 
observed for the shorter variant s-///Myb3 in the 
MBW combination. Both MMyb3 subvariants show a sig- 
nificant effect on both promoters in the MB combinations 
(Table 2, see also Figure 7) which is even more pro- 
nounced for Pc/zs4 (Table 2). HIMybl [19] appeared as the 
weakest activator of Pc/zs_Hl, although still an induction 
of 304 ± 47 pmol GUS/mg/min was observed for the 
MBW combination. Finally, although HlMyb7 shows 
some similarity to the other investigated Mybs (Figure 1), 
it is not an activator in any combination, neither for 
VchsJHl nor for Pc/zs4. 

HIMyb7 is a suppressor of promoter activation 

Although HIMybl is related to the other investigated 
Mybs, especially to A£Mybl2 and 7//Myb2, this hop- 
derived transcriptional regulator showed no ability to acti- 
vate the hop chs promoters (Figure 1). As we identified a 
transcriptional suppressor motif PDLNLELRIS, in the C 
terminal region of ///Myb7, the suppressive activity of 
HlMyb7 was investigated by assaying Pchs_Hl activity in 
transient expression systems containing MBW x HlMyb7 
combinations with ///Myb2, s-///Myb3 and A£Myb75 
(PAP1) (Figure 8). The addition of HlMyb7 was found to 
significantly decrease the reference GUS activity. A similar 



Table 2 Activation of promoters Pcfis_H1 and Pchs4 with 
selected Myb and hop TFs in various combinations 



Myb (M) 


AC 


M 


MB* 


MW** 


MB*W** 


Pd7S_H1 


H/Myb2 


FN646081 


(+) 


+ 


(+) 


++ 


s-H/Myb3 


AM501509 


0 


+ 


0 


++ 


/-H/Myb3 


AM501509 


0 


+ 


0 


0 


H/Myb1 # 


AJ876882 


0 


(+) 


n.d. 


(++) 


H/Myb7 


AB292244 


0 


0 


0 


0 


/\tMyb23 


NM_1 23397 


(+) 


(+) 


+ 


++ 


/4fMyb12 


NMJ30314 


++ 


++ 


++ 


0 


/\tMyb75 # 


NM_1 04541 


+ 


+ 


+ 


++ 


Pchs4 


H/Myb2 


FN646081 


0 


0 


0 


(+) 


s-H/Myb3 


AM501509 


+ 


+ 


+ 


0 


/-H/Myb3 


AM501509 


+ 


(++) 


0 


0 


H/Myb7 


AB292244 


0 


0 


0 


0 


/\tMyb23 


NM_1 23397 


0 


0 


(+) 


0 


/AfMyb12 


NM_1 30314 


+ 


0 


+ 


0 


/\tMyb75 $ 


NM_1 04541 


0 


(+) 


0 


(+) 



* B = H/DHLH2 (AC: FR751553) 
**W = H/WD4CM (AC: FN989721) 
$ M>AP1; 

# not tested for Pd?s4 

n.d. not determined 

(+) up to 10% of maximal level 

+ up to 70% of maximal level 

++ maximum activity observed; (++) low maximal activity; the maximal levels 
for M/promoter/combination/pmol GUS/mg/min: H/Myb2/Pd7sH1/MBW/5127; 
s-H/Myb3/Pd?sH1/MBW/1 192; I-H/Myb3/Pd?s4/MB/51 1; H/Myb1/Pd?s_H1/MBW/ 
304; AfMyb23/P_d7sH1/MBW/2036; /4fMyb12/P_d7sH1/MW/3108; /AfMyb75/ 
Pd7S_H1/MBW/5786 

0 = no increase of GUS activity over the background caused by M; no 
enhancement caused by MB or MW over M; no enhancement caused by MBW 
over M or over the binary combinations. 

suppression by HlMyb7 in the range of 50 to 90% was 
observed for binary combinations, as well as for the inde- 
pendent action of individual Myb genes on both chs pro- 
moters (data not shown), suggesting a more wide 
suppressive activity of this Myb. In order to analyze the 
possible role of the Myb-P boxes in this inhibition, the 
P2c/z5_Hl construct (Table 1) was used to analyze the influ- 
ence of HlMyb7 to the promoter activity after co-infiltra- 
tion with M2B2W1 x HMyb7, S-M3B2W1 x HlMyb7 and 
corresponding controls without HlMyb7. As the suppres- 
sive activity was also found to reach about 50% (data not 
shown), it is most likely that the Myb-P box does not play a 
major role. 

Influence of cloned TFs on the accumulation of 
anthocyanins in the heterologous system of petunia 
hybrida leaves 

The hop-derived TFs investigated in this study are all 
related to known transcription factors involved in 
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Figure 8 Suppressor action of H/Myb7 on Pd7s_H1 activation 

with TFs in various combinations. The combinatorial transient 

expression assay was performed and equal agrobacterium aliquots 

were mixed. For suppressor activity an aliquot of A tumefociens 

bearing 35S::HIMyb7 was added to the following combinations: M2, 

HIMyb2/HlbHLH2/HIWDR1; s-M3, s-HIMyb3/HlbHLH2/HIWDR1; PAP1, 

AtPAP 1/HlbHLH2/HIWDR 1 (white columns). Black columns, activity 

levels without HIMyb7. Instead of suppressor, LBA 4404 was added. 

C - control, infiltration with A tumefociens LBA 4404 without plant 

vector. Bars represent the confidence intervals at level a = 0.05. 
\ ) 

regulation of flavonoid pathways. Because of the signifi- 
cant activation of chs_Hl genes, the hop-derived TFs 
are expected to increase the pool of naringenin chal- 
cone, whereby they may stimulate downstream antho- 
cyanin biosynthesis and accumulation. In the lupulin 
glands, however, colored anthocyanins are not detectable 
as it is most likely that chalcone isomerase and further 
downstream enzymes of anthocyanin biosynthesis are not 
expressed. Instead, lupulin glands accumulate prenylated 
chalcones and other terpenophenolics composing the 
highly specific lupulin metabolome [52]. On the other 
hand, accumulation of naringenin chalcone in P. hybrida 
may lead to accumulation of anthocyanins. To further 
investigate the specificity of the cloned hop TFs in differ- 
ent combinations, as well as their potential to activate 
anthocyanin pathway, a transient TF expression system 
in P. hybrida was used. For comparison, AtVKVl was 
included as a regulator having the capability to co-activate 
the anthocyanin biosynthesis in petunia leaves providing 
co-expression of AtPAPl and chs_Hl genes arranged in 
tandem [16]). Hence, A£PAP1 or/ and hop chs_Hl activa- 
tors, e.g., s-///Myb3, ///Myb2, 7//bHLH2 and HIWDR1 
(Figure 7) in particular combinations were analyzed for 
possible anthocyanin induction. 

It was shown that only the co-infiltration of the chs_Hl 
gene (Figure 6B, construct III) with AtPAPl and HlbHLH2 
was found to strongly induce an accumulation of blue 
anthocyanins in the petunia leaves (Figure 9, for HPLC 



results see Additional file 5). This may be due to both an 
increased pool of chalcones and the activation of down- 
stream anthocyanin biosynthesis genes. Infiltration of 
either AtPAPl or HWHLH2 alone did not result in any 
detectable level of anthocyanins (Additional file 5). Based 
on retention times and UV-VIS spectra, the profile of the 
induced leaf anthocyanin metabolites was found to be 
qualitatively very similar to the profile of anthocyanins 
found in the blue corolla tissue of petunia flowers. How- 
ever, significant differences were found in quantitative 
composition. It is surprising that only minor amounts of 
anthocyanin metabolites were observed for the ternary 
combination A£PAP1B2W1 that led to dramatic activation 
of Pc/z5_Hl as shown above (Table 2). Neither the combi- 
nations with hop-derived Mybs, 7//Myb2 or s-///Myb3 that 
activated Pc/zs_Hl, i.e., 5-M3B2, S-M3B2W1 nor M2B2W1 
resulted in blue pigmentation in the petunia leaves. 
The metabolically highly active s-HlMyb3 [20], caused yel- 
low spots on the leaves in combination with 7//bHLH2 
(Figure 9). The hop TFs, ///Myb2, s-///Myb3 and WDR1 
were not able to co-induce anthocyanin biosynthesis in 
petunia leaves. 

Discussion 

Sequences and functional properties of lupulin gland- 
associated transcription factors and their potential role as 
regulators of the biosynthesis of the lupulin metabolome 

Four new transcriptional regulators, i.e., 7//Myb2, ///Myb7, 
7//bHLH2, and H/WDR1, were isolated and cloned from 
the Czech hop variety, Osvald's 72 in this study. All these 
TFs show a highly specific expression in the lupulin glands 
and display high similarity to established Myb, bHLH and 
WDR TFs, that are known regulators of the flavonoid bio- 
synthetic pathway in various plant species. Some of these 
regulators of flavonoid biosynthesis [24,61] form ternary 
MBW activation complexes like (PhAN2/ANl/ANll) in P. 
hybrida [29] and (AtTT2/TT8/TTGl) in A. thaliana [28] 
as well as recently described TF complexes in peas [30] and 
Lotus japonicus [31]. MBW complexes have been described 
to regulate not only pigmentation, but also the fate of epi- 
dermal cells, including the initiation and further develop- 
ment of trichomes [26,62]. Hop lupulin glands are 
glandular trichomes and, therefore, it was anticipated from 
the evolutionary point of view, that the combinatorial 
MBW machinery could be involved in the activation of 
genes related to the lupulin metabolome biosynthesis. 
Using a combinatorial transient expression system, it was 
shown that a M2B2W1 combination synergistically induced 
a very strong activation of native chsJM. This gene is cod- 
ing for a crucial hop "true" chalcone synthase endowed 
with the ability to produce naringenin chalcone with a high 
catalytical rate [16] and, therefore, it is likely playing an 
important role in the biosynthetic pathway leading to the 
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Figure 9 An example of anthocyanin pathway activation in 

Petunia hybrida leaves after infiltration with different TF 

combinations. The vector pLV-67 (see: Figure 6B, III) was infiltrated 

either in combination with s-HIMyb3 and HlbHLH2 (s-M3B2) or 

AtPAPl and HlbHLH2 (PAP1B2), mixed in equal amounts. The 

yellowish spots in the S-M3B2 variant and the anthocyanin spots in 

the PAP1B2 leaves developed ca. three days post-infiltration. The 

sites of infiltration are indicated by arrows, 
v . ) 



accumulation of prenylated chalcones in lupulin glands, 
including xanthohumol and desmethylxanthohumol (a 
direct precursor of the potent phytoestrogen 8-PN) [14,16]. 

The newly isolated HlMyb2 shows a typical sequence 
motif [DE]Lx2[RK]x3Lx6Lx3R within the R3 domain 
that is characteristic for a possible interaction with a 
member of the sub-group of Illf bHLHs [24,38,46]. This 
finding is in accordance with the ability of HlMyb2 to 
strongly activate Pc/zs_Hl in a M2B2W1 combination 
using a combinatorial transient expression system. On 
the other hand, clear synergistic effects were also observed 
for subvariants of the polyfunctional MMyb3 (cloned and 
characterized previously [20]), although the bHLH-inter- 
acting motif is missing in its amino acid sequence. Espe- 
cially the shorter subvariant s-///Myb3 resulted in a 
significant activation of the chs_Hl gene in both binary s- 
M3B2 and ternary S-M3B2W1 combinations. This sug- 
gests that 7//bHLH2 is able to interact with s-///Myb3 by 
either involving alternative hydrophobic residues in these 
positions exposed on the protein surface, or some alterna- 
tive way. According to our unpublished results, Electro- 
phoretic mobility shift assay (EMSA) showed some 
physical interaction between radioactively labeled VchsJHl 
and proteins extracted from leaves infiltrated with S-M3B2 
and S-M3B2W1 complexes. More experiments are neces- 
sary to clarify this interaction. The different length of the 
N-terminal domain of HlMyb3 has a crucial effect on its 
ability to activate Vchsjil in M3B2 and M3B2W1 combi- 
nations, as practically no synergism was observed for the 1- 
M3B2W1 combination. The divergent effects of overex- 
pressed 1- and s-///Myb3 TFs was observed earlier as 



reflected by changes in the composition of metabolites in 
petunia leaves, as well as in changes of the morphogenesis 
of A thaliana and P. hybrida transgenotes [20]. The func- 
tional analyses presented in Table 2 confirmed the com- 
plementation of hop lupulin gland-specific Mybs by 
selected A. thaliana Mybs in various binary or ternary 
combinations, suggesting not only sequential but also 
functional similarity of these TFs. For instance, AtVKVl in 
combination with MbHLH2 and 7//WDR1 resulted in a 
strong synergistic effect. Interplay between PAP1 and 
other components of MBW complexes in A. thaliana has 
been described previously [28]. In our experiments the 
highest synergism in the ternary MBW combination was 
also observed with A£Myb23, while the flavonoid regulator 
AtMybYl responded synergistically only to the binary BM 
and BW combinations, suggesting some specificity of pro- 
teimprotein interactions among these TFs. In addition to 
R2R3 Mybs, we confirmed the functional complementa- 
tion of HtWDRl by TTG1 from A. thaliana (not shown). 

The newly cloned HlMyb7 appears to be a potential 
negative R2R3Myb regulator of both chs _H1 and chs4 
genes in hop as it showed significant suppressor activity 
in a transient expression system. This is consistent with 
the identification of a PDLNLELRIS sequence in the C- 
terminal part of this protein conforming to the consen- 
sus pdLNLD/ELxiG/S sequence characteristic for sub- 
group 4 R2R3Mybs [63], which is conserved in 
repression region of the AtMyb4 repressor [42]. The 
repression caused by HlMyb7 was found to be indepen- 
dent on the Myb-P box positioned on the 5'region of 
Pc/zs_Hl, but it remains to be established whether it 
exerts its repressive activity via an interaction with 
MBW components, other proteins or binding to other 
elements in the promoter sequence. 

Upregulation of chs_Hl genes by both S-M3B2W1 and 
M2B2W1 combinations may increase the pool of available 
naringenin chalcone-precursors for further production of 
lupulin prenylchalcones co-determined by expression of 
downstream hop genes including prenyltransferase [64], 
and O- methyl transferase 1 [41]. While hop transformed 
by a single AtPAPl gene was characterized by a signifi- 
cantly higher production of anthocyanins, rutin and iso- 
quercetrin in cones and flowers in comparison to wild 
type plants [65], anthocyanins were not accumulated in 
the lupulin glands [52], most likely due to the lack of 
expression of the specific enzymatic machinery in these 
specialized structures. Conversely, hop TFs in S-M3B2W1 
and M2B2W1 combinations were unable to co-activate 
the anthocyanin pathway in P. hybrida as found for the 
A£PAP1B2 combination, thereby supporting the idea 
about functional specialization of lupulin-specific TFs. 
According to our unpublished results, neither the s- 
M3B2W1 nor the M2B2W1 combination was able to 
induce other promoters of hop genes like valerophenone 
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synthase or the O-methyl transferasel (Pomtl) [21], con- 
firming its specificity for the chs _H1 gene. 

Neither the ternary combinations S-M3B2W1, M2B2W1 
nor the Pc/zs_Hl -activating S-M3B2, M2B2 complexes 
were able to induce anthocyanin pigments in petunia 
leaves. The significantly elevated level of leaf metabolites 
and blue anthocyanins resulting from the A£PAP1B2 com- 
bination is probably due to an ability of A£PAP1B2 com- 
plex to activate also downstream genes in petunia, like 
chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), 
dihydroflavonol 4-reductase (DFR), and anthocyanidin 
synthase (ANS) [66]. Indeed, the metabolome and mRNA 
screenings shown that the overexpression of the A£PAP1 
gene is resulting in the activation of at least 25 petunia 
genes and also some new metabolites were produced [67] . 
Thus, further investigations are required to clarify the pro- 
file of hop genes that may be up-regulated by combinator- 
ial action of isolated hop TFs. 

Proposed regulation of hop chs genes as predicted from 
combinatorial analysis 

To analyze the functional properties of the newly cloned 
hop TFs it was opted to use a combinatorial transient 
expression system because hop transformation [68] and 
cultivation to reach flowering are very laborious and time- 
consuming especially for the investigation of the TFs in 
various combinations. While this transient expression 
assay is characterized by a negligible background signal i.e. 
signal originating from endogenous N. benihamiana TFs 
[21], it enabled us to efficiently compare the activation of 
different promoter constructs by various combinations of 
hop TFs. A comparison of truncated promoter constructs 
revealed a major dependence of Pchs_Hl activation by s- 
M3B2W1 and M2B2W1 combinations on a single Myb-P 
box characteristic for binding the PI TF of Zea mays [53]. 
However, other binding boxes such as the Myb-like box, 
the Myb-P-related box, the MYB IIG-like box and the 
G-binding sites characteristic for bHLH G-box binders 
driving the flavonoid pathway [23] also contribute to the 
overall activity. Therefore, each component of the promo- 
ter seems to impact on the degree of activation and it is 
the integral composition of the promoter that dictates the 
involvement of particular TFs in specific combinations. 
This is apparent on a comparison of two chs promoters 
both driving chalcone synthases having high expression 
in lupulin glands, i.e., chs_Hl and chs4 (see Figure 7, 
Table 2). A proposed schematic overview of the interac- 
tions of the investigated TFs with the different promoters 
is shown in Figure 10. The two promoters having a differ- 
ent spatial composition of similar potential Myb and 
bHLH binding boxes are activated differently by the same 
transcription factors. Pc/zs_Hl is predominantly activated 
not only by M2B2Wl/s-M3B2Wl ternary complexes, but 



also by B2W1 and MB binary combinations and by action 
of H/bHLH2 independent on hop B2 and Wl (Figure 10a, 
b,c,e). Single Myb and MW binary combinations did not 
result in activation of Pc/zs_Hl. This is in sharp contrast 
with Pc/zs4, where "direct activation" by HlMyb3 or the 
binary S-M3W1 combination is preferred (Figure 10d,e). 
Activation of Pc/zs4 by the M2B2W1 ternary complex is 
negligible (see Table 2), as well as any combination includ- 
ing 7//bHLH2 (Figure 10). Interestingly, both promoters 
can be negatively regulated by the HlMyb7 repressor 
(Figure 8; Figure 10). The general dependence of the com- 
binatorial action of TFs on the specific promoter composi- 
tion is illustrated in many studies summarized in recent 
reviews [23,24,69]. Although we did not observe any sig- 
nificant activation of Pc/zs_Hl with N. benihamiana TFs, 
we cannot exclude some influence of endogenous 
N. benthamiana TFs on the combinatorial regulation of 
Pc/zs_Hl in this heterologous system. Further study will be 
conducted to confirm proposed regulation pattern. The 
specific co-regulation of the lupulin metabolome biosynth- 
esis genes depends on the steady state levels of individual 
TFs forming the actual TF balance. It is, therefore, concei- 
vable that misbalancing of such complexes like M2B2W1 
by the downregulation of a specific TF complex compo- 
nent could lead to a significant decrease in gene expression 
and thus in the production of the lupulin metabolome. In 
support of this idea, it should be noticed that we recently 
detected a misbalanced expression of TFs in decolored 
petioles of diseased hop plants infected with hop stunt vir- 
oid (HSVd) (unpublished). 

Conclusions 

In the present study new TF homologues corresponding to 
the plant R2R3Myb, bHLH and WDR families were cloned 
using a cDNA library derived from hop (Humulus lupulus 
L.) lupulin gland tissue. The cloned ///Myb2 and -7, 
7//bHLH2 and H/WDR1 TFs display a high similarity to 
known TFs regulating the flavonoid biosynthesis pathway 
and were found to be highly specifically expressed in lupu- 
lin glands. The functional activity of these TFs was investi- 
gated using a combinatorial transient expression system in 
infiltrated leaf sectors of Nicotiana benthamiana. These 
experiments provided new insight into the complexity of 
the regulation of chs_Hl genes, as well as into the differ- 
ential activation of chs_Hl promoter and the lupulin-spe- 
cific chalcone synthase 4 promoter. The regulation 
involves variations of ternary, binary and independent TFs 
action. Complementation of hop and A. thaliana TFs was 
shown in particular combinations. The hop chs-activating 
TFs, H/Myb2 and 3, H/bHLH2 and H/WDR1 in various 
combinations did not enable the induction of the antho- 
cyanin pathway in P. hybrida suggesting the specialization 
of the MBW machinery in the specialized tissue of hop 
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Figure 10 Simplified scheme of Pd7s_H1 (I) and Pchs4 (II) activation with combinations of hop TFs as observed in heterologous 
system of N. benthamiana. The basic scheme was accommodated from Baudry et al. [27]. 3-D model of H/WDR1 in the scheme was obtained 
using SWISS-MODEL Workspace (see Methods), (a) the whole ternary complexes consisting of MYB, bHLH and WDR TFs; (b) bHLH and WDR 
binary combination; (c) MYB and bHLH combinations; (d) Myb and WDR combinations; (e) "direct" independent action of hop TFs that does not 
require complex with analyzed hop TFs; (f) inhibitory effect of HIMyb7. The figure is not in scale, TFs are assumed to interact with promoter 
binding sites. Binding boxes are not shown in the scheme. The weak signals close to background interactions are not considered in the scheme. 



glandular trichomes. HlMyb7 was characterized as an 
R2R3 repressor that could act as a potential co-regulator 
of the lupulin metabolome biosynthesis. 

Methods 

Plant material and plant cultivation conditions 

The Czech semi-early red-bine hop Osvald's clone 72 was 
grown under natural field conditions in the world hop col- 
lection of the Hop Research Institute in Zatec (Northern 
Bohemia). Petunia hybrida cv. Andrea and Nicotiana 
benthamiana plants were maintained in the green house 
at a temperature of 25 ± 3°C. Plants were grown under 
natural light conditions with supplementary illumination 
[170 umol m-2 s-1 PAR] to maintain a 16 h day period. 

Isolation and purification of lupulin glands, preparation 
of cDNA library and cloning of transcription factor cDNAs 

Lupulin glands were isolated from fully matured cones of 
the hop variety Osvald's clone 72 collected in the first half 
of August 2009 and 2010, two weeks before harvest using 
the procedure according to Nagel et al [41]. Isolated lupu- 
lin glands were used for RNA extraction, purification and 
hop cDNA library construction according to manufac- 
turers' instructions. Briefly, glandular tissue was frozen in 
liquid nitrogen and grinded in the Concert Plant RNA 
Reagent (Invitrogen). Then cell debris was removed by 
centrifugation and chloroform treatment, followed by 
RNA precipitation using isopropanol. Dissolved RNA pur- 
ification on columns with DNasel treatment was carried 
out according to the RNeasy Plant Mini Kit procedure 
(Qiagen), and subsequently mRNA was isolated with the 
Dynabeads mRNA Purification Kit (Invitrogen Dynal). 
Purified mRNA samples were used for cDNA library 



construction and lambda packaging, according to the 
ZAP-cDNA Synthesis Kit and the ZAP-cDNA Gigapack 
III Gold Packaging Kit instructions (Stratagene), 
respectively. 

The motifs for cloning of TFs cDNA were obtained by 
screening the Humulus lupulus sequences in both the 
GenBank EST (http://www.ncbi.nlm.nih.gov) and Tri- 
chOME (http://trichome.noble.org/trichomedb/) databases. 
In total, 258,472 ESTs were screened for R2R3Myb, bHLH 
and WDR analogues using the Advanced BLAST 2.0 script 
(http://blast.ncbi.nlm.nih.gov/). The obtained EST 
sequences were used for the design of primer sequences to 
be used for screening of the lupulin gland-specific cDNA 
library. High fidelity Pwo polymerase (Roche Molecular 
Biochemicals) was used for amplification. 

///Myb2 (clone 2962) was amplified from the library 
using primer combinations "Myb2 start" and "Myb2 stop", 
HIMybl (clone 2265) was amplified using primers "Myb7 
start" and "Myb7 stop" and for 7//bHLH2 (clone 3497) 
amplification we used the primers "bHLH start" and 
"bHLH stop". All primer sequences are listed in Additional 
file 6. For amplification of MWDR1(3140), a two-step 
amplification reaction was used. First, a mixture of cDNA 
fragments was obtained from the library using the primer 
combination "WDR start" and 3' oligo d(T)-anchor primer 
(573' RACE Kit 2 nd generation kit from Roche Diagnos- 
tics) and then a nested reaction was performed using a 
"WDRnes" primer together with "PCR WDR anchor" pri- 
mer. A fragment having the length of approximately 
1.6 kb was then cloned. 

cDNAs of A£Mybl2 and A£Myb23 were amplified 
from RNA of A. thaliana var. Columbia using the Titan 
One Tube RT-PCR (Molecular Biochemicals) including 
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a high fidelity Pwo polymerase (Roche Molecular 
Biochemicals) and primers derived from sequences 
described in the GenBank database (see Additional 
file 6). 

Genomic DNA preparation and genomic blots 

Hop genomic DNA was prepared by the CTAB method 
[70]. 5 ug DNA from each plant was digested with 50 
units of EcoRI, BamHl, Xbal, Xhol, Bgttl, or Pstl in 150 ul 
reaction mixes that were incubated overnight at 37°C. 
Digested DNA was subjected to electrophoresis and 
blotted as described previously [21] Pre-hybridization (2 h 
at 65°C in 30 ml) and hybridization (overnight at 65°C in 
20 ml) were carried out in buffer according to Church and 
Gilbert [71] using a hybridization oven. The membranes 
were hybridized to HIWDR1, HWHLH2 or HlMyb2 full 
length cDNA probes having a specific activity 5 x 10 7 
cpm/ug DNA. Hybridization signals were detected by 
phosphorimaging using a Typhoon 9200 imager (Amer- 
sham Pharmacia). 

RNA isolation and real-time PCR analyses 

For real-time quantitative PCR (RT qPCR) and for the 
reverse transcriptase-polymerase chain reaction (RT- 
PCR), total RNA was isolated from 100 mg of plant leaf 
tissue using Concert™Plant RNA Reagent (Invitrogen) 
following RNA purification and DNA cleavage on col- 
umns (RNeasy Plant Total RNA kit, Qiagen). The primer 
combinations for the RT qPCR of the hop TFs are listed 
in Additional file 6. All these primers had an annealing 
temperature close to 58°C and amplified unique, non- 
conserved regions of the TFs. Four micrograms of total 
RNA were reverse transcribed using oligo dT18 primer 
and Superscriptll reverse transcriptase (Invitrogen) at 
42°C for 60 min. A total of 5 |il of 50 x diluted cDNA 
was used for a 20 ul PCR reaction with 0.6 units of Hot 
Start Ex Taq polymerase (TaKaRa Bio), Taq buffer 1 x, 
dNTPs 200 uM each, SYBR Green 1:20,000 (Molecular 
Probes) and primers 375 nM each. All amplifications 
were carried out on a Bio-Rad IQ5 cycler for 40 cycles 
(94°C for 20 s, 59°C for 30 s, 72°C for 30 s) following an 
initial denaturation/Taq activation step (94°C for 5 min). 
The product size was confirmed by melting analysis and 
2% agarose gel electrophoresis. Data were analyzed and 
quantified with the Bio-Rad IQ5 Optical System version 
2.0 software. The abundance of a reference transcript, 
glyceraldehyde-3-phosphate dehydrogenase (GAPDH) 
[41,50,51], was estimated in parallel in each sample. 
GAPDH was amplified using the primer combination 
"HL-GAP-F1" and "HL-GAP-R1" as reported previously 
[21]. The relative values were standardized to GAPDH 
with the "Delta-delta method'and normalized to the sam- 
ple with highest expression (calibrator, set to 100%) 
according to Pfaffl [72] . The data points show the means 



± S.D. of two replicates of each PCR reaction. For RT- 
PCR amplifications we used cDNA primers as above and 
the high-fidelity Titan One Tube RT-PCR kit (Roche 
Molecular Biochemicals). 

Construction of plant expression vectors and preparation 
of chs_H1 promoter subvariants 

TFs HlMyb2, HlMyb7, HWHLH2 and HIWDR1 CDS 
were re- amplified with the high fidelity Pwo polymerase 
using the primer pair combinations (for the sequences 
see Additional file 6): "5'Myb2Apa" and "3 , Myb2Kpn"; 
"5'MYB7Xho" and "5'MYB7Xba; "bHLHKpn" and 
"bHLHBam"; "5'HlWDXho" and "3'HlWDXba". 

For the cloning of the A. thaliana TFs AtMybl2 and 
AtMyb23 we used the following combination of cloning 
primers: "3'Mybl2Xba" and "5 , Mybl2Xho"; "5'Myb23Xho" 
and "3'Myb23Xba". 

Prepared CDS fragments of HIWDR1, HlMyb7, 
AtMybl2 and AtMyb23 were treated with Xhol and Xbal; 
HlMyb2 with April and Kpnl and HlbHLH2 with Kpnl and 
BamHl. Gel-purified fragments were then fused to the 
35S CaMV promoter by cloning into the corresponding 
unique sites of the intermediary vector pLV68 [16]. The 
whole TF cassettes containing the 35S CaMV promoter 
were excised using restrictases Ascl and Pad and cloned 
in the plant vector pLV07 as described by Vrba and 
Matousek [73]. These vectors were introduced in the A. 
tumefaciens strain LBA 4404 by the freeze-thaw method 
[74] and used for leaf infiltration of N. benthamiana or P. 
hybrida leaves as described previously [16]. The ///Myb3- 
bearing plant vectors and A. tumefaciens strains were 
described previously [20]. 

The chsJHl promoter sequence was isolated earlier [14] 
and cloned to reference GUS vector previously [16]. The 
5'truncated variants (Table 2) were amplified using the 5' 
primers listed in Additional file 6: variant/primer: 1/ 
lPchsHlEco; 2/2PchsHlEc; 3/3PchsHlEco; 3/4PchsHl- 
Eco; 4/4PchsHlEco; 5/5PchsHlEco; 6/6PchsHlEco 
combined with 3' primer "PchslendXba". To prepare a 
mutated variant of the Myb-P box, the primer "mutP- 
myb5" was combined with the "PchslendXba" primer. 
The delta variant of VchsJHl was prepared using a combi- 
nation of "deltaPmyb5" and the "PchslendXba" primer. In 
this case variant 6 of the promoter served as a PCR tem- 
plate. The promoter Pchs4 (600 bp) was amplified from 
the chs4 gene [GenBank:AJ430353] [18] (clone 1746) 
using primers "Pchs4Eco" and "Pchs4Xba". The promoter 
regions were then fused to GFP/GUS genes by cloning in 
the unique EcoRI and Xbal restriction sites of the plant 
vector pBGF-0 [58] as described earlier [16]. The infiltra- 
tion with the vector pBGF-0 was used as negative control, 
while p35SBGF construct (clone 3771) containing 35S: 
GUS/GFP fusion cassette was used as p35S control. 
Summary of the strains of A. tumefaciens used and the 
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prepared T-DNA constructs can be found in the Addi- 
tional file 7. 

Combinatorial transient expression system and GUS 
assays 

For co-infiltration, agrobacteria suspensions were prepared 
at an OD 600 = 1.0 and equal volumes of each suspension 
were mixed prior to infiltration. N. benthamiana seedlings 
were sown in flats and grown under greenhouse condi- 
tions (22°C - 24° C) for 18 to 21 days. A tumefaciens cells 
carrying appropriate gene constructs were activated and 
injected into leaf blades as described previously [21]. The 
activity of ft -glucuronidase (GUS) was assayed in homoge- 
nized tissues of treated leaf sectors 3-4 d following the 
infiltration using a fluorometric assay [75] modified by 
Matousek et al. [21]. The amount of the released fluores- 
cent dye 4-methylumbelliferone (MU) was measured using 
the VersaFluor™ fluorometer (Bio-Rad) with excitation at 
365 nm and emission measured at 455 nm. The fluorom- 
eter was calibrated with a fresh preparation of MU (100 
nM) as standard. GUS activity was expressed in pmol MU 
.mg 1 fresh tissue .min 1 . At least three independent 
repeated experiments were performed in each experimen- 
tal variant. The data were statistically analyzed using 
Microsoft Excel 2010 software. 

Analysis of anthocyanin metabolites 

Leaves of P. hybrida were lyophilized prior to analysis of 
anthocyanins. After grounding up lyophilized plant 
material in a mortar, samples of 10-20 mg were 
extracted in eppendorfs using 1 mL of methanol/acetic 
acid/water (9/1/10, v/v/v) and, after sonication, the 
extraction mixtures were kept for 16 h at 4°C. After 
centrifugation at 18000 rpm for 10 min, the supernatant 
was analyzed by high performance liquid chromatogra- 
phy (HPLC) using a Waters Alliance 2695 Separation 
Module (Waters, Zellik, Belgium) equipped with a 996 
photodiode array detector (PDA). The column was an 
Agilent Zorbax Bonus RP (Agilent, Diegem, Belgium, 
3.5 |im, 150 x 3.0 mm), while the injection volume was 
10 |iL. A linear gradient elution at a flow rate of 0.7 
mL/min was applied using increasing ratios of solvent B 
(acetonitrile/acetic acid/phosphoric acid/water, 25/20/ 
1.5/53.5 v/v/v/v) in solvent A (water/phosphoric acid, 
98.5/1.5, v/v) according to the following schedule: 0 
min: 25% B in A, 12 min: 35% B in A, 18 min: 48% B in 
A, 37 min: 100% B in A. After 5 min of washing at 
100% of solvent B, the column was re-equilibrated at 
25% B in A for 10 min. Chromatograms at 530 nm were 
extracted from the 3D-data. Peak integrations were car- 
ried out using standard parameters. Other secondary 
metabolites in the leaves of P. hybrida were analyzed 
according as previously described [21]. 



Bioinformatic analyses 

For sequence comparisons and cluster analyses we used, 
in addition to cloned TFs H/Myb2, H/Myb7, HlbHUYl 
and ///WDR1, Myb, bHLH and WDR amino acid 
sequences from nucleotide and protein databases as pro- 
vided by NCBI services (see Additional file 2, Additional 
file 3 and Additional file 4). Sequence analyses were car- 
ried out with Phylogeny Pipeline (http://phylogeny. 
lirmm.fr/phylo_cgi/simple_phylogeny.cgi) [76]. The trees 
were generated using the PhyML v. 3.0 (http://www.atgc- 
montpellier.fr/phyml) [77]. The generated data were 
visualised by FigTree v.1.3.1. (http://tree.bio.ed.ac.uk/ 
software/figtree). Some sequence comparisons were car- 
ried out with DNASIS for Windows, version 2.6 (Hitachi 
Software Engineering Company, Tokyo, Japan). Amino 
acid classes described by [78] were used for protein com- 
parisons. Plant C/5-acting Regulatory DNA elements 
(PLACE) database [79] and Promoterlnspector option of 
the Genomatix (Eldorado Gene2Promoter) software 
(http://www.genomatix.de) were used to screen for TFs 
binding sites. A comparative protein analysis of the 3-D 
structures of the NR2R3 regions of H/Myb2 and H/Myb3 
was performed by means of the SWISS-MODEL Work- 
space (http://swissmodel.expasy.Org//SWISS-MODEL. 
htm) [80] using the loading template lh88C.pdb. Align- 
ments of 3-D structures and structural analyses were per- 
formed using Swiss-Pdb Viewer v3.7b2 [81]. 

Additional material 



Additional file 1: List of R2R3Myb TFs included in the phylogenetic 
tree presented in Figure 1 and alignment of amino acid sequences 
within R2R3 domain. 

Additional file 2: Comparative analysis of the 3-D protein structures 
of the N-terminal/R2/R3 regions of H/Myb2, I-H/Myb3 and s-H/Myb3 

The theoretical structures were calculated and portrayed against the 
template 1h88C.pdb using the SWISS-MODEL Workspace. The alignments 
of the 3-D structures were performed using the Swiss-PdbViewer v3.7b2 
(see Methods). Structures are presented as single ribbons, positions of R2 
and R3 repeats are shown on alignment s-/-//Myb3/l-/-//Myb3; the 
positions of essential residues forming the bHLH-binding site of H/Myb2 
are shown on the structure alignment for s-H/Myb3/H/Myb2. The 
positions of hydrophobic amino acids substituting the bHLH-binding 
residues in s-H/Myb3 are between brackets. The positions of the largest 
structural deviations are indicated on the structures by filled and hollow 
arrows. N and C are N- and C- terminus, respectively. 

Additional file 3: List of bHLH TFs included in the phylogenetic tree 
presented in Figure 2A and alignment of amino acid sequences 
within bHLH domain. 

Additional file 4: List of WDR TFs included in the phylogenetic tree 
presented in Figure 3Aand alignment of amino acid sequences. 

Additional file 5: Relative levels of anthocyanin pigments in petunia 
leaves (%) infiltrated with different hop TFs or >4fPAP1 compared to 
anthocyanins in petunia dark blue corollas performed by HPLC 
analysis. 

Additional file 6: Oligonucleotide primers used in this study. Table 
of oligonucleotide primers used in this study, arranged according the 
purpose for which they were designated. 
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Additional file 7: List of Agrobacterium tumefaciens strains and 
vectors. Table of A tumefaciens strains and vectors used in the analysis 
of lupulin gland transcription factors from R2R3Myb, bHLH and WDR 
families. 

Additional file 8: Histochemical and quantitative analysis of GUS 
activity driven by modified variants of Pchs_H1 by hop TFs 
(supplement to Table 1). 
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8-PN: 8-prenylnaringenin; ANS: anthocyanidin synthase; B2: /-//bHLH 2; DFR: 
dihydroflavonol 4-reductase; EMSA: electrophoretic mobility shift assay; F3H: 
flavanone 3-hydroxylase; GADPH: glyceraldehyde-3-phosphate 
dehydrogenase; GUS: B-glucuronidase; HPLC: high performance liquid 
chromatography; HSVd: hop stunt viroid; CHI: chalcone isomerase; CHS: 
chalcone synthase; M2: H/Myb2; M3: H/Myb3; OMT: 0-methyl transferasel; P: 
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qPCR: real-time quantitative PCR; TF: transcription factor; VPS: valerophenone 
synthase; W1: H/WDR1; X: Xanthohumol. 
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